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ABSTRACT 



Context. FU Orionis-type objects (FUors) are low-mass young eruptive stars that probably represent an evolutionary phase charac- 
terized by episodic periods of increased accretion rate from the circumstellar disk to the star. Theory predicts that a circumstellar 
envelope, the source of continuous mass infall onto the disk, is necessary for triggering such accretion bursts. 

Aims. We intend to study the spatial and velocity structure of circumstellar envelopes around FUors by means of molecular line 
observations at millimeter wavelengths. We target three prototypical FUors as well as an object possibly in a pre-outburst state. 
Methods. We present archival PdBI interferometric observations of the J=1~0 line of I3 CO at 110.2 GHz. For three of our targets, 
these represent the first millimeter interferometric observations. The data allow the study of the molecular environment of the objects 
on a spatial resolution of a thousand AU and a velocity resolution of 0.2 km s -1 . 

Results. Strong, narrow 13 CO(1-0) line emission is detected from all sources. The emission is spatially resolved in all cases, with 
deconvolved sizes of a few thousand AUs. For V1057 Cyg and V1331 Cyg, the emitting area is rather compact, suggesting that the 
origin of the emission is a circumstellar envelope surrounding the central star. For V1735 Cyg, the 13 CO emission is offset from the 
stellar position, indicating that the source of this emission may be a small foreground cloud, also responsible for the high reddening 
of the central star. The 13 CO emission towards V1515Cyg is the most extended in the sample, and apparently coincides with the 
ring-like optical reflection nebula associated with V1515 Cyg. 

Conclusions. We suggest that millimeter interferometric observations are indispensable for a complete understanding of the circum- 
stellar environment of FUors. Any theory of the FUor phenomenon that interprets the geometry of the circumstellar structure and its 
evolution using single beam measurements must be checked and compared to interferometric observations in the future. 
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1. Introduction 

FU Orionis-type objects, or shortly FUors, constitute a 
small group of young stars characterized by large out- 
bursts in visible light , attrib uted to highly enhanced accretion 
dHartmann & Kenvoni [l996). During these outbursts, accretion 
rates from the circumstellar disk to the star are in the order of 
10~ 4 M Q /yr, three orders of magnitude higher than in quiescence. 
Enhanced accretion is often accompanied by enhanced mass 
loss: most FUors have optical jets, molecular outflows, and op- 
tically visible ring-like structures that are sometimes explained 
by expanding shells thrown off during previous outbursts. 

During a single, century-long outburst, as much as 0.01 M 
of material can be dumped onto the stellar surface. Thus, the in- 
ner disk needs to be replenished after each outburst, possibly by 
material from the infalling envelope. Recent theoretical studies 
show that the continuous infal l from the envelope is al so nec- 
essary to trigger the outbursts (Vorobvov & BasuL [20061) . After 
many outbursts, the envelope vanishes, and the object finally en- 
ters a state of permanently low accretion. This general paradigm 
of the evolution o f young, low-mass stars was invoked by 
Oua nz et al.l (120071) to explain the observed diversity of FUors: 
while some objects are still deeply embedded (e.g. L1551 IRS 5, 
V1735Cyg), others have already cleared away part of their en- 
velopes (e.g. V1057Cyg, V1515Cyg). Statistics show that prob- 
ably all low-mass young stars undergo FUor-like phases during 



their evolution, implying that FUors might be the clue objects to 
study envelope evolution and dispersal. 

Interferometric observations of molecular line emission have 
been successfully used to probe small-scale structure of molec- 
ular material around young stellar objects. For example, for 
a sam ple of low-mass YSOs in T a urus dHogerheiide et"aTl 
119981). Serpens dHogerheiide et al.L 119991) . and Ophiuchus 



dvan Kempen et aUl2009l) . it was found that most YSOs are sur- 
rounded by compact envelopes of up to a few thousand AU in 
radius, well traced by 13 CO and C ls O. Condensations, inho- 
mogeneities in the envelopes can be seen in HCO + and 13 CO. 
The HCO + and HCN molecules trace the walls of the outflow, 
while SiO and SO emission originates from shocked material 
in the outflow. These observations provided detailed kinemati- 
cal picture of the envelopes (for molecu lar studies of some in- 
dividual young stellar objects, see e.g. lJ0rgensen et al.l [2004; 
iMatthews et alJl2006l;lBrinch et al.ll2009l) . 

Despite being "accretion laboratories", there are relatively 
few molecular gas observations published for FUors. Single-dish 
observations often show the presence of strong CO emission to- 
wards FUors, and line pro files sometimes indicate the presence 
of molecular outflows (e.g. lHartmann & Kenvonlll996L and ref- 
erences therein). However, due to the large single dish beam, it is 
not possible to study the spatial distribution of the emission us- 
ing these data. Out of the 9 FUors listed in Hart mann & Kenvoni 
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(1996), millimeter interfe rometric observation s are only avail- 
able for L1551 IRS 5 (e.g.|Momoseetall[i99l)- 

In this paper, we present interferometric observations of the 
molecular emission from four FUor-type objects in Cygnus: 
V1057Cyg, V1331Cyg, V1515Cvg , and V1735Cvg. With the 
excep tion of V1331Cyg dLevreaulll [19881: iMcMuldroch et al.L 
1 1993h . these data represent the first millimeter interferometric 
observations for our targets. We use the 13 CO(1-0) line emis- 
sion to perform a high spatial and spectral resolution study of the 
gaseous material around FUors. We analyze the spatial and kine- 
matic structure of the 13 CO gas, and check whether the emission 
can be associated with the circumstellar envelopes. Our results 
can be compared with those obtained for normal YSOs. Such a 
comparison may also contribute to the on-going debate whether 
all young stars undergo an eruptive phase in their early evolution 
or FUors are atypical objects. 

2. Observations 

We reduced unpublished 13 CO observations of our targets ob- 
tained with the Plateau de Bure Interferometer (PdBI) on May 
30 and 31, 1993 (program ID: C057, PI: D. Fiebig). The ob- 
servations were carried out in snapshot mode (about 1 hour 
on-source correlation time per object), using four antennas 
(4D1 configuration), and the baselines ranged from 24 to 64 m. 
The receiver was tuned to the I3 CO(1-0) line at 110.2GHz 
(lower sideband), the channel spacing was 78 kHz (0.21 km/s). 
At this wavelength, the single dish HPBW is 45". Bright 
quasars (3C454.3, 2005+403, 2037+51 1, 2200+420, 0648-165, 
and 0727-115) were observed to enable RF bandpass, phase, 
and amplitude calibration. The phase center was 20 h 58 m 53 s .7 
+44°15'27'/9 for V1057Cyg, 21 h 01 m 09 s .3 +50°21'42'.'2 for 
V1331Cyg, 20 h 23 m 47.9 +42°12'24'.'8 for V1515Cyg, and 
21 h 47 m 20 s 7 +47 °32'04'.'l for V1735Cyg. 

Data reduction was done in the standard way with CLIC, a 
GILDAS-based application written especially for the reduction 
of PdBI dat^]. The rms phase noise was less than 18° for all 
observations. We estimate a flux calibration accuracy of «15%. 
MAPPING was used to create naturally weighted dirty maps 
in an area of 64" x 64", which were then deconvolved using 
the Clark CLEAN method. The synthesized beam is approxi- 
mately 7" x 6" while the rms noise in the channel maps is ^0. 15 
Jy/beam. For the exact phase centers and beam parameters, see 
the log of observations in Table Q] 



3. Results 

The 13 CO(1-0) spectra of our sources, integrated over the 
whole 64"x64" maps, are plotted in Fig. Q] Strong 13 CO(l- 
0) emission was detected towards all sources with high signal- 
to-noise ratio. In all cases, we found an emission line coincid- 
ing with the syst emic velocity of each source (from optica l and 
near-IR spectra. iHerbid 1 1 9771 iKenvon & Hartmannlll989l and 
IChavarriaet all 19791 determined heliocentric radial velocities of 
-14 ±2, -15 ±2, -12 ±2, and -10 ±2kms~ 1 , coresponding to 
LSR velocities of 1.9, 0.6, 5.2, and 3.9 kms 1 for V1057Cyg, 
V1331 Cyg, V1515Cyg, and V1735Cyg, respectively). We col- 
lapsed a few channels centered on the peak emission, and cal- 
culated visibility amplitudes as a function of uv radius (Fig. [2]), 
and produced velocity-integrated maps (Fig.0. The LSR veloc- 
ities where the emission peaks are listed in Table Q] Velocity- 
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Fig. 1. I3 CO(1-0) spectra of our targets obtained with the PdBI. 
Fluxes were integrated over an area of 64" x 64" centered on the 
phase center. 



integrated line fluxes calculated for the emission areas visible in 
Fig. [3] are also listed in Table Q] 



3. 1 . Line profiles 

The line profile of V1057Cyg is the broadest in our sam- 
ple (FWHM»;2.1kms _1 ), and is centered at 4.6kms~ 1 . This 
is consistent with the velocity seen in single dish 1 2 CO(l- 
0), | 2 CO(2-l) and 13 CO(1-0) data by iBechis&Lol (fT97l 
and iLevreaull d!988l) . V13 31Cyg shows a na r row 1 3 CO(1-0) 
line (FWHM=s0.7kms _1 ). IMcMuldroch et all (1 19931) presents 
Owens Valley interferometric observations of VI 331 Cyg in 
the same line, as well as in other transitions and other iso- 
topes of CO. The velocities of these lines are consistent with 
ours. V1515Cyg shows two separate, narrow emission com- 
ponents (FWHMa;0.6kms _1 ) at 5.1kms _1 and at 11.9kms~'. 
Both of these components are also visible i n sing le dish 12 CO(3- 
2) and 13 CO(2-l) data from lEvans etal.1 (QJ94). The spectrum 
of V1735Cyg shows a narrow, sing le-peaked 13 CO(1-0) line 
(FWHMw0.7kms _1 ). The shape and position of this line co- 
incides well with that of the 13 C O(2-l) line detected in single 
dish data bv lEvans et all (11994 . The 12 CO(2-l) line from the 
same paper shows self-absorption at this velocity, and broad line 
wings indicating outflow activity. We plotted position-velocity 
diagrams along different angles through our sources, but except 
for V 1 5 1 5 Cyg (which will be discussed later), we found no sig- 
nificant velocity gradients. 



http://www.iram.fr/IRAMFR/GILDAS 



3.2. Spatial extent 

In Fig. [2] we plotted the visibility amplitudes as a function of 
uv radius. Since these are snapshot observations, the Mv-space is 
only partially sampled; there are ranges of uv radii where no vis- 
ibility information is available. However, even with this coarse 
sampling, the graphs show a decline of amplitude with uv ra- 
dius, indicating that our sources are spatially resolved. Simple 
2D Gaussian fits to the maps presented in Fig. [3] also indicate 
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Table 1. Summary of PdBI observations. Distances are from S andell & Weintraubl (120011) . The FWHM of the fitted Gaussians are 



deconvolved sizes. 
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Fig. 2. Visibility amplitudes as a function of uv radius. Data were 
binned in 1 m- wide bins and the error bars indicate the dispersion 
of data points within one bin. 
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Fig. 3. 13 CO(1-0) maps of our targets, integrated in the velocity 
ranges indicated above the images. White asterisks indicate op- 
tical positions from SIMBAD. Beam sizes are ^6x7". The noise 
level is cr=0.10Jykms _I ; the solid contours are 2<r, Act, 6<x, 
the dashed contours are 0, -2<x, -4cr. 



that our targets are spatially resolved. The deconvolved sizes and 
PAs of the Gaussian fits are listed in Table [TJ Gaussian models 
fitted to the visibilities give sizes within +2" and PAs within 
+ 10° to those determined from the Gaussian fits to the images. 
Thus, we consider these numbers as representative uncertainties 
of the values given in Table Q] We emphasize that these fits are 
not physical envelope models, and thus the obtained parameters 
should be regarded only as rough quantitative estimates of the 
size of the emitting region. 

V1057Cyg and V1331 Cyg are the most compact sources in 
our sample, with FWHM of about 10" x5". Moreover, in both 
cases, the emission in centered precisely on the optical position 
of the stars. V15 15 Cyg shows the most extended emission in our 
sample. The stellar position seems to be located at the southern 
tip of a bright, slightly curved, elongated filament, but extended 
emission can be seen both to the south and to the north of the 
star (»20" across). A map produced for the channels around the 
emission line at 11.9kms _1 reveals that this separate velocity 
component comes from a compact area ^13" to the northwest 
(P.A. 330°) of the optical source. The bulk of the emission in the 
VI 735 Cyg region comes from an area with a FWHM of «9" 
centered 2'.'9 to the northwest with respect to the optical position 
of the star. Moreover, there is also some tentative extended emis- 
sion towards the northwest about 14" from the stellar position. 



3.3. Gas masses 



Following lScoville et al ] (119861) . the total H? mass can be calcu- 
lated from the observed I3 CO(1-0) line fluxes by 



M Hl = 2.39 x 10" 



(r x + 0.89) T.3 CO 



^kpc 



-5.31/7* 



1 - e- T X( l3 CO) 



f 



S v dv M e 



where Tx is the excitation temperature, r is the optical depth, 
D is the distance of the source, X( 13 CO) is the molecular 
abundance relative to H2, and J S v dv is the velocity-integrated 
line flux in units of Jykms 1 . Using 7^= 50 K (as determined 
from dust continuum observations by ISandell & Wein traub 
2001), assuming the lines to be optically t hin (t «:1), and 
taking X( 13 CO)=1.6xl0~ 6 (lLanger & PenziasL 1 19931) . the total 
gas masses are 0.12, 0.06, 0.43, and O.1OM for V1057Cyg, 
V1331Cyg, V1515Cyg, and V1735Cyg, respectively. Of 
course these numbers should be considered lower limits if the 
lines are actually optically thick, or if part of the 13 CO emis- 
sion is resolved out by the interferometer. However, these val- 
ues are in good agreement with total masses derived from 
sin gle-dish 850 um dust continuum maps of comparable size 
bv ISandell & W eintraub: dlOOlh . which were 0.10, 0.13, 0.15, 
and 0.42M o for V1057Cyg, V1331Cyg, V1515Cyg, and 
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V1735Cyg, respectively. This suggests that our interferometric 
observations probably recover most of the 13 CO emission. 

3.4. 2.7mm continuum 

Although the main focus of our data analysis is the study of the 
13 CO emission, the PdBI observations could also be used to look 
for 2.7 mm continuum emission by excluding channels around 
the lines and collapsing the remaining channels. The resulting 
continuum maps have a typical rms noise of 2 mJy/beam. We 
detected 2.7 mm continuum emission only for one of our tar- 
gets, VI 331 C yg. This is not surpris i ng be cause, based on the 
observations of lSandell & Weintraubl (1200 ll) . out of our four tar- 
gets, VI 331 Cyg is the brightest also at 1.3 mm and at 850 jum. 
The 2.7 mm emission we detected is not resolved, its position 
is consistent with the optical stellar position, and the measured 
total flux is 12+2 mJy. Assuming optically thin emission, dus t 
opacity of *i.3 mm =0.01 cm 2 g _1 (Ossenkop f & Henningl 1 19941) . 
an emissivity law of k~A~ 1 and dust temperature of 50 K 
dSandell & Weintrau b. 2001), this flux corresponds to a total 
(gas+dust) mass of O.19M , similar to what we obtained from 
the line flux of V1331Cyg. Using 6mJy as a 3cr upper limit 
for the 2.7 mm continuum flux of the other sources, we can 
give upper limits of 0.11, 0.32, and 0.26M G for V1057Cyg, 
V 15 15 Cyg, and VI 735 Cyg, respectively. 

4. Discussion 

In the following subsections, we discuss what our interferomet- 
ric 13 CO(1-0) data can add to what we already know about the 
circumstellar environment of our targets based on optical im- 
ages, sub-millimeter continuum maps, CO line observations, and 
SED analysis from the literature. We will also discuss whether 
the detected 13 CO emission can be associated with the en- 
velopes. 

4.1. V1 057 Cyg 

After its outburst in 1969-70, an eccent ric ring-like n ebulosity of 
1' x l'.5 appeared around V1057Cyg dHerbigUl977l) . The mon- 
itoring of this ring in the subsequent years indicated that it is 
fading together with the central star without any changes in the 
structure. This implies that the ring is a reflection nebula: a pre- 
existing structure illuminated by the outbursting star, and not 
material ejected from the central source by the outburst. From 
the line wings of sing le dish 12 CO data, iLevreaultl (fl988h and 
lEvans et all {1994) concluded the presence of a molecular out- 
flow on a similar spatial scale. The sur roundings of VI 057 Cyg 
was mapped by Sandell & Weintraub (2001), whose 850yum 
continuum image indicates the presence of a rather compact but 
resolved («4") source coinciding with the star and a fainter, 
north-south oriented filament. 

In our 13 CO(1-0) data we do not detect emission either from 
the ring-like reflection nebula of from the north-south filament. 
However, our interferometric observations definitely resolve the 
central source: the deconvolved Gaussian FWHM at the distance 
of V1057Cyg corresponds to about 5800x3100AU. The lack 
of observable line wings, the relatively narrow line profile, and 
the compact spatial appearance of the emission points to quies- 
cent gas, probably situated in the outer parts of a circumstellar 
envelope. 

The infrared SED of the source (e.g. Abraham e t all 120041) 
suggests a continuous temperature distribution in the system. 



This supports our conclusion that there is an envelope associ- 
ated with and heated by the central star. The presence and mea- 
sured size of the envelope is consistent with th e model predic- 
tion of 7000 AU proposed by Gre en et all (12006). They also sug- 
gested that the envelope has a rather large conical cavity which 
explains the modest far-infrared excess. The relatively low ex- 
tinction along the line-of-sight and the detectability of the 10 fim 
silicate emission feature implies a close to pole-on geometry for 
the source, in accordance with its relatively symmetric shape 
in our interferometric observations. The circumstellar mass of 
O.12M derived from ou r 13 CO observations, consist ently with 
the dust mass estimate bv lSandell & Weintraubl d2001l) . is rather 
high, significantly exceeding the typical disk masses of T Tauri- 
type stars. 

4.2. V1331Cyg 

V 1331 Cyg is not a FUor in its present state, but - based on 
the similarity of its spectrum to that of VI 057 Cyg prior to 
maximum light - is probably in a pre-o utburst state, or be- 
tween outbursts (McMul droch et afl. Il993l) . Si ngle-dish and in- 
terf erometric 12 CO and 13 CO observations by iLevreau lil (fl98l 
and IrvlcMuldr och et all d 19931) revealed a complex circumstel- 
lar environment containing a molecular outflow approximately 
along the line of sight, a flattened gaseous envelope of about 
6000 x 4400 AU in size, and a gaseous expanding ring of 
about 41 000x28 000 AU. The latter coincides well with the 
large r ing-shaped optical reflection nebula seen bv lOuanz et all 
d2003). 

Our interferometric 13 CO observations of VI 331 Cyg, con- 
sistently with the s imilar angular resolution measurements of 
iMcMuldroch et all (Il993l) . reveal a rather compact structure to- 
wards the star with a deconvolved Gaussian FWHM of about 
5300 x 3600 AU at the distance of V1331Cyg. This core is 
also visible in continuum emission in the 850 /mi maps of 
ISandell & Weintraubl (gOOl), although it is only marginally re- 
solved (<6" or 3300 AU). We suggest that both the 13 CO gas 
emission and the dust thermal emission is originated in a cir- 
cumstellar envel ope, probably the same fl attened structure that 
was proposed by McMuldroch et al. ( 1993) to explain the nature 
of the 13 CO emission. It is probable that the dust in this flattened 
envelope scatters the stellar li ght and gives rise to the inner ring 
observed in optical images by Ouanz et al 1 d2007l) . 

The presence of a circumstellar envelope around VI 331 Cyg 
is also consistent with the SED, which exhibits significant in- 
frared excess indicating substantial amounts of circumstellar 
material (e.g. [A braham et all 12004]). The geom e try pr oposed by 
McMuldr och et alj (1 19931) and by lOuanz et all (2007) suggests 
that we see the inner part of the system through a conical cav- 
ity filled with a pole-on molecular outflow. However, we detect 
I3 CO emission neither from the outflow nor from the outer, large 
expandin g ring, in agreement wit h the interferometric observa- 
tions of McMuldroc h et al.1 dl993l) . 

4.3. V1515Cyg 

VI 551 Cyg brightened in optical light slowly during the 1940s 
and 1950s (Herbig, Il977l) . Photographic plates from this time 
show a bright, narrow arc of nebulosity extending to the north 
and west of the star. Later images from the 1970s to the present 
day still show this northern arc, but a brighter nebulosity is also 
visible to the south and west, together forming a nearly com- 
plete circular ring with a diameter of «T6" (corresponding to 
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Right Ascension 

Fig. 4. The optical and millimeter environment of V1515Cyg. 
The grayscale image in the background is an SDSS r-band im- 
age; the gray contours are the same 13 CO(1-0) contours as in 
Fig. [3] White asterisk indicate the optical stellar position. 



16000 AU, see also the Sloan Digital Sky Survey r-band image 
from 2003 in Fig. 0). The fact that the size of this ring did not 
change significantly for 60 years but the brightness distribution 
did, suggests that it is a reflection nebula similar in nature to that 
around V1057Cyg. 

Our 13 CO(1-0) measurement shows arc-shaped emission, 
which, plotted over the optical image in Fig. [4] clearly coin- 
cides with the ring-shaped reflection nebula. The similar mor- 
phology of the 13 CO emission and the scattered light, as well as 
the matching radial velocities of the 13 CO emission and the star 
suggest that the detected molecular gas is physically associated 
with the FUor. S imilarly to VI 057 Cyg, the 12 CO line profile 
of V1515Cyg in lEvans et al.1 (11994 indicate the presence of a 
molecular outflow. Our 13 CO(1-0) data, however, show no high- 
velocity line wings; the emission line is narrow and thus proba- 
bly traces quiescent, not outflowing material. We detect a small 
(<0.5kms~') velocity difference between the northern and the 
southern part of the ring, the southern part being slightly more 
red-shifted. However, the ring is clearly not expanding. 

A map of the 850/vm dust continuum from 

Sand ell & Weintraubl (1200 ll) shows faint extended emis- 
sion towards V1515Cyg and also at about 10" to the north 
and northwest to the star. Although these observations were 
made with a relatively large (15") beam, the results are not 
inconsistent with a potentially arc-shaped dust emission. It is 
thus possible that the ring of material that is responsible for the 
optical reflection nebula and for the 13 CO emission also emits 
in dust continuum. 

Although the observed molecular emission does not exactly 
peak towards the star, some dust and gas must be located in the 
vicinity of V 15 15 Cyg forming a circumstellar envelope. This 
claim is suggested by the SED of the FUor, which was mod- 
eled with a d isk+e nvelope geometry bv lTurner et al.l (1 19971) and 
iGreen etaD J2006). However, this envelope does not show up 
as a separate localized peak in our 13 CO map. We measured the 
highest gas mass towards this object within our sample, which 



can probably be explained by the contribution from both the cir- 
cumstellar envelope and from the arc-shaped feature. Thus the 
0.43 M is an upper limit for the envelope mass. 

4.4. V1 735 Cyg 

V 1735 Cyg h ad its outburst some time in the 1950s or the 1960s 
(lEiiasLll978l) . The star is associated with faint patches and fila- 
ments of reflecte d light within 1'. S ubmillimeter dust continuum 
observations by ISandell & Weintraubl (1200 lh indicate the pres- 
ence of two sources: one associated with (albeit slightly offset 
from the optical position of) VI 735 Cyg, and a brighter one lo- 
cated about 20-24" to the northeast (the deeply embedded Class I 
protosta r V1735 Cyg SMn. Single-dish 12 CO and 13 CO obser- 
vations dLevreaultl 1 19831; [Richardson et all 119851; lEvans et all 
IT991 revealed a complex molecular gas structure in the vicin- 
ity of these sources. The 12 CO emission is extended in the ar- 
cminu te spatial scale, an d broad line wings indicate outflow ac- 
tivity. lEvans et all (119941) concluded that both VI 735 Cyg and 
SMI drive molecular outflows. 

Similarly to the other three FUors in our sample, our 13 CO 
interferometric observations of V1735Cyg show narrow emis- 
sion with no line wings. This again points to quiescent gas in 
some kind of envelope. Our map shows a compact peak close 
to VI 735 Cyg, but with a definite offse t of 2"9 to the n orthwest. 
Based on multi-filter optical images, iGoodrichl (1 19871) discov- 
ered that V1735 Cyg itself is much more reddened than the sur- 
rounding reflection nebulosity. He speculates that the reason for 
this may be a dark cloud in front of V1735Cyg. In this sce- 
nario, the dark cloud may be offset from the star, as long as it 
causes enough interstellar extinction to explain the observed col- 
ors. Based on the fact that the 13 CO emission we detect is offset 
from the stellar position, we propose that a significant part of the 
emission is coming from this foregro und dark cloud. Th is sug- 
gestion is supported by the results of lOuanz et al . (2007), who, 
based on the analysis of mid-infrared ice features, concluded 
that the extinction for VI 735 Cyg might be caused by ices some- 
where in the line of sight to the source, rather than material re- 
lated to the young star. The foreground structure, however, may 
not be completely unrelated to the star because of their identical 
radial velocities. 

For continuum ob servation s such as those presente d in 
IWeintraub etail d!991l) and in ISandell & Weintraubl d2001l) . the 
dominant source is not the FUor, but VI 735 Cyg SMI. This 
source also emits molecular line emission, as evidence d by th e 
single-dish 12 CO(3-2) map presented by Eva ns et alj (Q994). 
Interestingly, we did not detect V1735CygSMl in 13 CO. The 
reason for the non-detection may be partly due to the fact that 
V1735CygSMl is separated from V1735Cyg by about 20", 
thus it is at the very edge of our primary beam. Additionally, 
the submillimeter source is very extended, thus the interferome- 
ter may filter out most of its emission. 

5. Summary and Conclusions 

In this paper, we present interferometric observations of the 
13 CO(1-0) line of four well known FUors, young stellar objects 
characterized by large optical outbursts due to enhanced disk ac- 
cretion. For V1057 Cyg, V1515Cyg, and V1735Cyg, these rep- 
resent the first millimeter interferometric data published so far. 
This makes it possible to study the gas distribution on a few thou- 
sand AU spatial scale. 

Although all of our sources are known to drive molecu- 
lar outflows, as evidenced by the high-velocity line wings of 
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12 CO, our data suggest that the 13 CO(1-0) emission traces 
quiescent gas. With the exception of V1515Cyg, the size of 
the emitting region is within a few thousand AU, consistent 
with typical circumstellar envelope sizes. Gas masses calculated 
from our 13 CO line fluxes a nd from dust continuum data from 
ISandell & Weintraubl J200l are also consistent. This indicates 
that the 13 CO emission seen towards our sources is mostly origi- 
nated from a relatively compact circumstellar envelope, or, as in 
the case of V1735 Cyg, possibly from a small foreground cloud. 

All of our sources are surrounded by reflection nebulosities, 
which are probably pre-existent structures illuminated by the 
brightened central source. With the exception of V 15 15 Cyg, we 
detect no 13 CO emission from these structures. For V1515 Cyg, 
the 13 CO emission coincides with the ring-shaped optical reflec- 
tion nebula. This indicates that the source is surrounded by a ring 
of material that on one hand scatters the optical light of the cen- 
tral star, and on the other hand emits at millimeter wavelengths. 
An important consequence is that attributing unresolved, single- 
dish millimeter fluxes of V 15 15 Cyg to the circumstellar enve- 
lope may not be entirely correct. 

Based on the appearance of the lOyum silicate feature, 
Ouanzetal 1 (f2007h defined two categories of FUors. They argue 
that objects showing the feature in absorption (e.g. V1735Cyg) 
are younger, still embedded in a circumstellar envelope. 
Objects showing the silicate band in emission (e.g. V1057Cyg, 
V1331 Cyg, and V1515 Cyg) are more evolved, with direct view 
on the surface layer of the accretion disk. However, in case 
of V1735Cyg, the millimeter emission and the absorption at 
shorter wavelengths may be due to a foreground cloud, and 
consequently the object may be more evolved (less embed- 
ded) than it appears. Our 13 CO data, especially the case of 
V 15 15 Cyg and VI 735 Cyg, demonstrate that millimeter emis- 
sion and consequently absorption at other wavelengths may not 
be necessarily or exclusively associated with circumstellar en- 
velopes. Thus, large beam, single-dish data alone are probably 
not enough to obtain a complete picture of the circumstellar en- 
vironment of FUors. We suggest that any theory of the FUor phe- 
nomenon that interprets the geometry of the circumstellar mate- 
rial and its evolution using single beam measurements must be 
checked and compared to interferometric observations in the fu- 
ture. Millimeter observations with facilities such as ALMA for 
both the dust continuum and for the line emission of molecules 
such as 13 CO, C ls O, HCO + , HCN, etc. will be a fruitful direction 
in future studies of FUors. 
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